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The bound state generating functional is constructed in gauge theories. This con- 
struction is based on the Dirac Hamiltonian approach to gauge theories, the Poincare 
group classification of fields and their nonlocal bound states, and the Markov- Yukawa 
constraint of irreducibility. The generating functional contains additional anomalous 
creations of pseudoscalar bound states: para-positronium in QED and mesons in 
QCD in the two gamma processes of the type of 7 + 7 = ttq + para-positronium. 
The functional allows us to establish physically clear and transparent relations be- 
tween the perturbative QCD to its nonperturbative low energy model by means of 
normal ordering and the quark and gluon condensates. In the limit of small current 
quark masses, the Gell-Mann-Oakes-Renner relation is derived from the Schwinger- 
Dyson (SD) and Bethe-Salpeter (BS) equations. The constituent quark masses can 
be calculated from a self-consistent non-linear equation. 

Is dedicated to the 60-th anniversary of the birth of Professor S.I. Vinitsky 

1. INTRODUCTION 

At the beginning of the sixties of the twentieth century Feynman found that the naive 
generalization of his method of construction of QED fails in the non-Abelian theories. The 
unitary S-matrix in the non-Abelian theory was obtained in the form of the FP path integral 
by the brilliant application of the theory of connections in vector bundle [1]. Many physicists 
are of opinion that the FP path integral is the highest level of quantum description of the 
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gauge constrained systems. Anyway, the FP integral at least allows us to prove the both 
renormalizability of the unified theory of electroweak interactions and asymptotic freedom 
of the non-Abelian theory. However, the generalization of the FP path integral to the bound 
states in the non-Abelian theories still remains a serious and challenging problem. 

The bound states in gauge theories are usually considered in the framework of representa- 
tions of the homogeneous Lorentz group and the FP functional in one of the Lorentz-invariant 
gauges. In particular, the "Lorentz gauge formulation" was discussed in the review [2] with 
almost 400 papers before 1992 on this subject being cited. Presently, the situation is not 
changed, because the gauge-invariance of the FP path integral is proved only for the scatter- 
ing processes of elementary particles on their mass shells in the framework of the "Lorentz 
gauge formulation" [3]. 

In this paper, we suggest a systematic scheme of the bound state generalization of FP- 
functional and S-matrix elements. The scheme is based on irreducible representations of the 
nonhomogeneous Poincare group in concordance with the first QED description of bound 
states [4, 5]. This approach includes the following elements. 

i) The concept of the in- and out- state rays [6] as the products of the Poincare represen- 
tations of the Markov- Yukawa bound states [7-10]. 

a) The split of the potential components from the radiation ones in a rest frame. 

in) Construction of the bound state functional in the presence of the radiation compo- 
nents. This functional contains the triangle axial anomalies with additional time derivatives 
of the radiation components. 

iv) The joint Hamiltonian approach to the sum of the both standard time derivatives and 
triangle anomaly derivatives. 

All these elements together leads to new anomalous processes in the strong magnetic 
fields. One of them is the two gamma para-positronium creation accompanied by the pion 
creation of the type of 7 + 7 = 7ro+ para-positronium. 

Within the bound state generalization of the FP integral, we establish physically clear 
and transparent relations between the parton QCD model and the Numbu-Jona-Lasinio ones 
[11-13]. Below we show that it can be done by means of the gluon and quark condensates, 
introduced via the normal ordering. 

The paper is organized as follows. Sec. 2 regards the Poincare classification of in- and out- 
states. In Sec. 3, the Dirac method of gauge invariant separation of potential and radiation 
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variables is considered within QED. Sec. 4 is devoted to the bound state generahzation of 
the Faddeev-Popov generating functional. In Sec. 5, we discuss the bound state functional in 
the presence of the radiation components, which contains the triangle axial anomalies with 
additional time derivatives of these radiation components. Sec. 6 is devoted to the axial 
anomalies in the NJL model inspired by QCD. In Appendix A, the bound state functional in 
the ladder approximation is considered. In Appendix B, the BS equations are written down 
explicitly and discussed. 

2. BOGOLIUBOV-LOGUNOV-TODOROV RAYS AS IN-, OUT-STATES 

According to the general principles of quantum field theory (QFT), physical states of the 
lowest order of perturbation theory are completely covered by local fields as particle-like 
representations of the Poincare group of transformations of four-dimensional space-time. 

The existence of each elementary particle is associated with a quantum fields tp. These 
fields are operators defined in all space-time and acting on states \V, s) in the Hilbert space 
with positively defined scalar product. The states correspond to the wave functions \E'ci(x) = 
{0\i/ja{x)\V , s) of free particles. 

Its algebra is formed by generators of the four translations = id^ and six rotations 
M^'^ = i[x^j_du — Xyd^^. The unitary and irreducible representations are eigen-states of the 
Casimir operators of mass and spin, given by 

V^\V,s) = ml (1) 
-wl\V,s) = s{s + l)\V,s), (2) 
Wp = ^ex^upV'M^". (3) 

The unitary irreducible Poincare representations describe wave-like dynamical local exci- 
tations of two transverse photons 

^5) -) = / (0 E 71=^^(^)" [e^(— ) A,-,, + e-(— ) A," J . (4) 

Two independent polarizations ei^b)a are perpendicular to the wave vector and to each other, 
and the photon dispersion is given by = v^l?. 

The creation and annihilation operators of photon obey the commutation relations 
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The bound states of elementary particles (fermions) are associated with bilocal quantum 
fields formed by the instantaneous potentials (see [7-9]) 

M(x,y) ^ M(z\X) (5) 



H 



where V ■ X = ujhXq — PX, = {ujh, P) is the total momentum components on the mass 



shell (that is, uh = y Mjj + V ), and 

X = ^—, z = x-y. (6) 

are the the total coordinate and the relative one, respectively. The functions F belongs to the 
complete set of orthonormalized solutions of the BS equation [4] in a specific gauge theory, 
0^(1^, q~^) are coefficients treated in quantum theory as the creation (+) and annihilation 
operators (see Appendix B). 

The irreducibility constraint, called Markov- Yukawa constraint, is imposed on the class 
of instantaneous bound states 

z^V,M{z\X) ^ ^z^^J^Miz\X) = 0. (7) 

In Ref. [6] the in- and out- asymptotical states are the "rays" defined as a product of these 
irreducible representations of the Poincare group 

(out| = |, \m) = \l[Vj,sj). (8) 

J J 

This means that all particles (elementary and composite) are far enough from each other to 
neglect their interactions in the in-, out- states. All their asymptotical states (out| and |in) 
including the bound states are considered as the irreducible representations of the Poincare 
group. 

These irreducible representations form a complete set of states, and the reference frames 



Ml 



are distinguished by the eigenvalues of the appropriate time operator 

i,\V,s) = '^\Vj,s), (9) 
where the Bogoliubov-Logunov-Todorov rays (9) can include bound states. 
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3. SYMMETRY OF S-MATRIX 

The S-matrix elements are defined as the evolution operator expectation values between 
in- and out- states 

Min^ = (pu^ , (10) 

P—inv,G—inv P— variant P— variant, G—inv P— variant 

where the abbreviation "G-inv", or "gauge-invariant", assumes the invariance of S-matrix with 
respect to the gauge transformations. 

The Dirac approach to gauge-invariant S-matrix was formulated at the rest frame £° = 
(1,0,0,0) [14-16]. Then the problem arises how to construct a gauge-invariant S-matrix in 
an arbitrary frame of reference. It was Heisenberg and Pauli's question to von Neumann: 
"How to generalize the Dirac Hamiltonian approach to QED of 1927 [14] to any frame?" 
[10, 15, 16, 27]. The von Neumann reply was to go back to the initial Lorentz-invariant 
formulation and to choose the comoving frame 

£j = (1,0,0,0) = £^, £/^ = £.£ = l (11) 

and to repeat the gauge-invariant Dirac scheme in this frame. 

Dirac Hamiltonian approach to QED of 1927 was based on the constraint-shell action [14] 

(12) 



QED_q' 



where the component is defined by the scalar product Aq = A ■ i oi vector field and 
the unit time-like vector i^. 

The gauge was established by Dirac as the first integral of the Gauss constraint 

dt^-^ = 0, t={x.i). (13) 
OAq 

In this case, the S-matrix elements (10) are relativistic invariant and independent of the 
frame reference provided the condition (9) is fulfilled [9, 18]. 

Therefore, such relativistic bound states can be successfully included in the relativistic 
covariant unitary perturbation theory [18]. They satisfy the Markov- Yukawa constraint (7). 

This framework yields the observed spectrum of bound states in QED [5], which corre- 
sponds to the instantaneous potential interaction and paves a way for constructing a bound 
state generating functional. The functional construction is based on the Poincare group 
representations (51) with £^ being the eigenvalue of the total momentum operator of instan- 
taneous bound states. 



4. QED 

4.1. Split of potential part from radiation one 

Let us formulate the statement of the bound state problem in the terms of the gauge- 
invariant variables using QED. It is given by the action [16] 

W[A^p,i>] = Jd'x [-^(F^,)2 + ^(^)^(A)-mO)^] , (14) 

V^(A) = d^-icA^, yj = v^-Y , 

Ff^u = d^A,-d,A^ . (15) 

Dirac defined these gauge-invariant variables by the transformations 

E 4^^ = ^^H = v\A\{A,^i-d^{v{A\r\ (16) 

a=l,2 ^ 

= v{AM,, (17) 

where the gauge factor is given by 

v\J\\ = exp{ze J dt'aoit')} (18) 

ao[A] = ^mA{t,x) . (19) 

Here the inverse Laplace operator acts on arbitrary function /(t, x) as 

i/(,x) S -I j d»,ii^^ (20) 

with the kernel being the Coulomb potential. 
Using the gauge transformations 

a^ = ao + 9oA =^ = exp[ieA(to, x)]t;[A] exp[-ieA(t, x)] , (21) 

we can find that initial data of the gauge-invariant Dirac variables (16) are degenerated with 
respect to the stationary gauge transformations 

[A^] = Af[A] + 9,A(to,x), = exp[2eA(to,x)]^^[A,^] . (22) 

The Dirac variables (16) as the functionals of the initial fields satisfy the Gauss law constraint 

5o(9,Af(t,x)) = . (23) 
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Thus, explicit resolving the Gauss law allows us to remove two degrees of freedom and to 
reduce the gauge group into the subgroup of the stationary gauge transformations (22). 

We can fix a stationary phase A{to, x) = $o(x) by an additional constraint in the form of 
the time integral of the Gauss law constraint (23) with zero initial data 

diA^ = ^ A$o(x) = 0. (24) 

Dirac constructed the unconstrained system, equivalent to the initial theory (14) 

W* = W\sw/5Ao=o[A^ = a:, = r] (25) 
'{A*f-Bf , 1 .,1 



+ t; Jo IT Jo - jT A* + ri^d - m)r 



2 2"'' A 



where 



M = ^0^^"' (26) 



a=l,2 



= E^jkdjAl (27) 

are the electric and magnetic fields, respectively. 

In three-dimensional QED, there is a subtle difference between the model (25) and the 
initial gauge theory (14). This is the origin of the current conservation law. In the initial 
constrained system (14), the current conservation law 9ojo = diji follows from the equa- 
tions for the gauge fields, whereas a similar law 8^]^ = dij* in the equivalent unconstrained 
system (25) follows only from the classical equations for the fermion fields. This difference 
becomes essential in quantum theory. In the second case, we cannot use the current conser- 
vation law if the quantum fermions are off mass-shell, in particular, in a bound state. What 
do we observe in an atom? The bare fermions, or dressed ones (16)? Dirac supposed [14] 
that we can observe only gauge invariant quantities of the type of the dressed fields. 



4.2. Bilocal fields in the Ladder Approximation 

The constraint-shell QED allows us to construct the relativistic covariant perturbation 
theory with respect to radiation corrections [19]. Recall that our solution of the problem 
of relativistic invariance of the nonlocal objects is the choice of the time vector 
operator with eigenvalues proportional to the total momenta of bound states [20, 21]. In this 
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case, the relativistic covariant unitary S-matrix can be defined as the Feynman path integral 

J*] = (*| J DiP*Dij*e'^^^^*'^*^-^'^'\*) , (28) 

where 



+ 11 A(V'(z/)^(x))/CW(z^ I Xmxmy))], (29) 

and the symbol 

(*|...|*) = / Y[DA*/^Sm _ . . (30) 
^ j 

stands for the averaging over transverse photons. Here by definition p = d^'j^, and /C^^-* is 
the kernel 

JC^'\z^ I X) = ^V{z^)6{z ■ £)/, (31) 

where z and X are the relative and total coordinates (6). The potential V{z-^) depends only 
on the transverse component of the relative coordinate with respect to the time axis i. The 
requirement for the choice of the time axis (9) in bifocal dynamics is equivalent to Markov 
- Yukawa condition (7). 

Apparently, the most straightforward way for constructing a theory of bound states is the 
redefinition of action (29) in terms of the bifocal fields by means of the Legendre transfor- 
mation [22] 

^ / d^xd%{ilj{y)ip{x))}C{x,y){ilj{x)il!{y)) = 

= _i J d^xd^yM{x,y)K:-\x,y)Mix,y) + (32) 
+ / d'^xd^y{ij{x)^{y)),M{x,y) 

where IC~^ is the inverse kernel /C given by Eq. (31). Following Ref. [23], we introduce the 
short-hand notation 

d^xd^ytlj{y)i!{x){ip - le/i* - m°)6^^\x - y) = {^j^j, -G^^) , (33) 
dVy{ip{x)ilj{y))M{x,y) = {'^ip,M). (34) 
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After quantization over fermion fields (here A^^ is the number of colors equal to 3), the 
functional (28) takes the form 

Z?[s,r,J*] = (*| jY[DMe'^^ff^^^+'^^ff^^^\*), (35) 

where 

W,fj[M] = tr [log(-G^i + A^)] - ^iM,JC-'M) (36) 
is the effective action, and 

S,ff[M] = is*s*,iG~/-Mr') (37) 

is the source term. The effective action can be decomposed as 

1 °° 1 

WefflM] = --{M, JC^'M) + z Yl (38) 

n=l 

Here $ = Ga-M, "J?'^, etc. mean the following expressions 

$(a;, y) = GaM = / d^zGA{x, z)M{z, y), 



= j d^xd^y^{x,y)^{y,x), (39) 

$^ = y d'^xd'^yd'^z^{x,y)^{y, z)^{z,x), etc. 

As a result of such quantization, only the contributions with inner fermionic lines (but not 
the scattering and dissociation channel contributions) are included in the effective action, 
since we are interested only in the bound states constructed as unitary representations of 
the Poincare group. 

4.3. The anomalous creation of Para-Positronium in QED 

The effective bound state functional in the presence of radiation fields contains a triangle 
anomaly decay of positronium rjp with an additional time derivative of these fields 

W,ff = W{A*) + W{7]p) 



W{A*) = j d^x[Gp7]pAB, 

W, = Jd^xl^llvl-Mlvp'-id^VP?]], 
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where Bi are the magnetic field component (27), and the parameter of the effective action 
is given by 



2a ftsJ^^\ v^«'/' «V2 
Gp = 



where a = 1/137 is the QED couphng constant, and 

Fp 



(40) 



(41) 



is the positronium analogy of the pion weak-coupling constant F^^ discussed in the next 
Section. 

The product CpAiBi is obtained from the triangle diagram shown in Fig. 1. 



rXAAAAAAA/ 



V 




Figure 1. The standard triangle diagram, used for calculating Cp, the parameter of the effective action. 

The Hamiltonian of this system is the sum of the Hamiltonians of the free electro- magnetic 
fields and the positronium ones r]p{x) and the interaction 



Weff = J dtd^x\^E,Ai + Pr,f]p-n} 



n,.t = Clr^pE^B, + ^^B} . (42) 

The anomalous processes of creation of the positronium pairs in the external magnetic field 
at the photon energy value E^i ~ mg (see Fig. 2) are described by the differential cross section 



da_ ^ a'E', _ /2meY ^ ^a''E'^ _ /2mey 
dn 7rl28F|\/ V / 512m4 y V / 
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Figure 2. The new diagram for the anomalous processes of creation of the positronium pairs in the 
external magnetic field that follows from the Hamiltonian anomaly, given by Hint in Eqs. (42). 

4.4. The Schwinger QEDi+i 

The Schwinger two dimensional QEDi+i was considered in the framework of the Dirac 
approach to gauge theories distinguished by the constraint-shell action [24]. 

This constraint-shell action has an additional time derivative term of the gauge field that 
goes from the fermion propagator in the axial anomaly. This anomalous time derivative term 
changes the initial Hamiltonian structure of the gauge field action 

W^Scwinger = / dtdx \ ^f]l + CsVsA + ^ i (44) 



dtdx 



Psfjs + eA-^ + CsVsE - Cf-^ 



(45) 



e 



Cs = ^. (46) 

Finally, an additional Abelian anomaly given by the last term in Eq. (45) enables us to 
determine the mass of the pseudoscalar bound state [24]. In (5-E-D(i+i), it is the well-known 
mass of the Schwinger bound state 

AM^ = -. 

TT 

The Schwinger model justifies including of the similar additional terms in the 4-dimensional 
QED. 
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5. NON-ABELIAN DIRAC HAMILTONIAN DYNAMICS IN AN ARBITRARY 

FRAME OF REFERENCE 



In order to demonstrate the Lorentz-invariant version of the Dirac method [14] given by 
Eq.(12) in a non-Abehan theory, we consider the simplest example of the Lorentz-invariant 
formulation of the naive path integral without any ghost fields and FP-determinant 



HdA 



(47) 



We use standard the QCD action t/^, ■?/'] and the source terms 



W 



d X 



F,\ = doAt - doAld + g^AlAl ^ At - VfA^ 
S = Jd^x [A^ + # + ^r/] , = ig^^. 



(48) 
(49) 
(50) 



There are a lot of drawbacks of this path integral from the point of view of the Faddeev- 
Popov functional [1]. They are the following: 

1. The time component Aq has indefinite metric. 

2. The integral (47) contained the infinite gauge factor. 

3. The bound state spectrum contains tachyons. 

4. The analytical properties of field propagators are gauge dependent. 

5. Operator foundation is absent [25]. 

6. Low-energy region does not separate from the high-energy one. 

All these defects can be removed by the integration over the indefinite metric time com- 
ponent Af^^^ = A- i, where is an arbitrary unit time-like vector: = 1. If £° = (1, 0, 0, 0) 
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then Afj^i'^ = ^o- In this case 



n 



dtVf{A*)Af = 0, 



a\2 



,4 

d X — — 



d'^xip 



hodo-1 j{dj + A*)- m 



d^xdSj^ix) 



;5^{x - y) 



ab 



(51) 

(52) 

(53) 
(54) 
(55) 

(56) 



The infinite factor is removed by the gauge fixing (52) treated as an antiderivative function 
of the Gauss constraint. A*"" denotes fields A" under gauge fixing condition (52). It becomes 
homogeneous \/i^{A*)A*^ = because Aq is determined by the interactions of currents (56). 
It is just the non-Abelian generalization [10, 26-28] of the Dirac approach to QED [14]. In 
the case of QCD there is a possibility to include the nonzero condensate of transverse gluons 

The Lorentz-invariant bound state matrix elements can be obtained, if we choose the 
time-axis i of Dirac Hamiltonian dynamics as the operator acting in the complete set of 
bound states (9) and given by Eqs.(6) and (7). This means the von Neumann substitution 
(11) given in [15] 

Z[f]^Z[i]^Z[i] (57) 
instead of the Lorentz-gauge formulation [1]. 



6. AXIAL ANOMALIES IN THE NJL MODEL INSPIRED BY QCD 

6.1. Formulation of the NJL model inspired by QCD 

Instantaneous QCD interactions are described by the non-Abelian generalization of the 
Dirac gauge in QED 

1 



5'inst = J d^xq{x) {ip - rh^)q{x) - \ j d'^xd^yJoix) 



L(V,(A*)) 



-8 (x-y) 



Joiy) (58) 
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where j^ix) = (j{x) —'yoq{x) is the 4-th component of the quark current, with the Gell-Mann 
color matrices A" (see the notations in Appendix A). The symbol m° = diag(m°, mj], m°) 
denotes the bare quark mass matrix. 

The normal ordering of the transverse gluons in the nonlinear action (56) V^^A^V^'^Aq 
leads to the condensate of gluons 

g^f-^^f-^'^{ArAf*) = 2g^[Nl - l]<5''^<5.,Cgi,on = M^S'^,, (59) 

where 

{AfAf) = 2Cgi,on%5"'. (60) 

This condensate yields the squired effective gluon mass in the squared covariant derivative 
ydft^ftydc^c ^. ^dh^br^dc^c^ . j^M'^A'^A'^ of coustraiut-shell action (56) given in Appendix 
A. The constant 

^ gluon 



(27r)32Vk2 

is finite after substraction of the infinite volume contribution, and its value is determined 
by the hadron size like the Casimir vacuum energy [29]. Finally, in the lowest order of 
perturbation theory, this gluon condensation yields the effective Yukawa potential in the 
colorless meson sector 

and the NJL type model with the effective gluon mass M^. While deriving the last equation, 
we use the relation 



a=l 



colorless 



in the colorless meson sector. 

Below we consider the potential model (58) in the form 

^mst = j d^xqix) {%p - m°)g(x) - ^ ^ d^xd%3l{x)V{x^ - y^)b{{x - y) ■ i)j^{y) (62) 

with the choice of the time axis as the eigenvalues of the bound state total momentum, in 
the framework of the ladder approximation given in Appendix A. 
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6.2. Schwinger-Dyson equation: the fermion spectrum 

The equation of stationarity (A6) can be rewritten from the SD equation 

S(x — y) = 'm^5^^\x — y) + ilC{x, y)GY.{x — y) . (63) 

It describes the spectrum of Dirac particles in bound states. In the momentum space with 
S(/c) = j (i^xS(x)e*'^'^' for the Coulomb type kernel, we obtain the following equation for the 
mass operator ( S ) 

m) = ^' + lJ - ^^)^^s(^)^' (64) 

where Gj]{q) = — '^{q))"^ ■, Y_{k^) is the Fourier representation of the potential, /c^ = 
kfj. — ifj.{k ■ i) is the relative transverse momentum. The quantity S depends only on the 
transverse momentum S(/c) = because of the instantaneous form of the potential 

V_{k-^). We can put 

S,(g) = ^a(q) cos2t;a(q) = Ma(q). (65) 

Here Ma(q) is the constituent quark mass and 

o t \ ^a(q) 
cos2t;a(q) = , ^ (66) 

determines the Foldy - Wouthuysen type matrix 

S'a(q) = exp[(q7/g)t;a(g)] = cost;a(g) + (qj/q) smv^{q) (67) 

with the vector of Dirac matrices 7 = (71,72,73) and some angle Va{q). The fermion spec- 
trum can be obtained by solving the SD equation (64). It can integrated over the longitudinal 
momentum qo = {q ■ i) in the reference frame = (1, 0, 0, 0), where q-^ = (0, q). By using 
Eq. (67), the quark Green function can be presented in the form 



= [lo^-E.i<l^)S-\q^)]-' = 

(68) 



go - E^{q^) + ie qo + E^{q^) + ie 



where 



AS2)a(?^) = ^a(g^)AS2)(0)5a"'(9^), A|2)(0) = (1 ± /)/2 (69) 
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are the operators separating the states with positive (+-E'a) and negative {—Eg) energies. As 
a result, we obtain the following equations for the one-particles energy E and the angle v 
with the potential given by Eq. (61) 

E,{k^) cos2v,{k^) = ml + ^j ^V{k^ - q^) cos2v,{q^). (70) 

In the rest frame £^ = (1, 0, 0, 0) this equation takes the form 

M.ik) = ml + ^J -0^Vik - q) cos2t;,(g). (71) 

By using the integral over the solid angle 



M2 + (k-q)2 J ^ M2 + P + g2 _ 2A;ge kq + {k - qf 

and the definition of the QCD coupling constant as = 4:7rg'^, it can be rewritten as 

oo 

M(k)-m' + ^ f dq^MM=ln^^l±^^±^ (72) 
Mg[k) - + y d^^ Mliq) + q^ + (fc - g)^" ^^^^ 

The suggested scheme allows us to consider the SD equation (71) in the limit when the 
bare current mass equals to zero. Then the ultraviolet divergence is absent, and, hence, 
the renormalization procedure can be successfully avoided. 

This kind of nonlinear integral equations was considered in the paper [30] numerically. 

The solutions show us that in the region q <^ Mg the function cos2t;a is almost constant 

cos2t;a — 1, whereas in the region q ^ Mg the function cos2t;a(g) decays in accordance with 

the power law {Mg/qY~^^ . The parameter /3 is a solution of the equation 

cot(/37r/2) 3 

= 2' (73) 

lying in the range < /3 < 2. This equation has two roots for < as < S/n, the first, 
belonging to the interval < /3i < 1, and the second, related to the first one by = 2 — 
At as = 3 /n, the two solutions merges into /3 = 1, and there is no root for larger values of 
the coupling constant. Equation (73) can be obtained by means of linearization of Eq. (71) 
within the range q ^ Mg, because in this range Ma_{q) <^ q. Thus, the solution for cos2fa(g) 
is a reminiscent of the step function. This result justifies the estimation of the quark and 
meson spectra in the separable approximation [21] in agreement with the experimental data. 
Currently, numerical solutions of the nonlinear equation (72) are under way, and the details 
of computations will be published elsewhere. 
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6.3. Spontaneous chiral symmetry breaking 

As discussed in the previous section, the SD equation (71) can be rewritten in the form 
(72). Once we know the solution of Eq. (72) for Ma(g), we can determine the Foldy- 
Wouthuysen angles t^a, (a = u,d) for u-,d- quarks with the help of relation (66). Then 
the BS equations in the form (BIO) 

M^L-(p) = [i?„(p) + i?rf(p)]Lnp)-y"^Z(p-q)i^^(q)[c-(p)c-(g)+es-(p)s-(g)], (74) 

M^LUp) = [Euip) + E,{p)]Ll{p)-j -^nP-^)Limc^{p)c^{q)Hs'^{p)s^{q)] (75) 

yield the pion mass M^^ and wave functions i^i(p) and L2(p). Here m^jm^ are the current 
quark masses, E;^ = a/p^ + M^{p), (a = u,d) are the u-,d- quark energies, ^ = {p<l)/pq, and 
we use the notations 

E{p) = E,{p) + i?b(p) , (76) 
c^(p) = cos[t;a(p) ± Vb{p)] , (77) 
s^{p) = sm[v^{p) ±Vb{p)] . (78) 

The model is simplified in some limiting cases. Once the quark masses m„ and are 
small and approximately equal, then Eqs. (71) and (74) take the form 

m. = A4(p) - i y -^np - q) COS 2«„(5), (79) 

= B.(p)Lr(p)-i/|f3Hp-q)iJ(q). (80) 

Solutions of equations of this type are considered in the numerous papers [31-35] (see also 
review [30]) for different potentials. One of the main results of these papers was the pure 
quantum effect of spontaneous chiral symmetry breaking. In this case, the instantaneous 
interaction leads to rearrangement of the perturbation series and strongly changes the spec- 
trum of elementary excitations and bound states in contrast to the naive perturbation theory. 

In the limit of massless quarks m„ = the left-hand side of Eq. (79) is equal to zero. The 
nonzero solution of Eq. (79) implies that there exists a mode with zero pion mass M^- = in 
accordance with the Goldstone theorem. This means that the BS equation (80), being the 
equation for the wave function of the Goldstone pion, coincides with the the SD equation 
(79) in the case of = M^^ = 0. Comparing the equations yields 

rn(^ Mujp) cos2vu{p) 
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where the constant of the proportionahty F^^ in Eq. (81) is called the weak decay constant. 
In the more general case of massive quark ^ M^^ ^ 0, this constant is determined from 
the normalization condition (B17) 

with = 3. In this case the wave function Ll{p) is proportional to the Fourier component 
of the quark condensate 

n=Nc 



It, ^YC p 

Cquark = J] (^^ ^)^n , Y) ) = 4iVc / 
n.=l 



Using Eqs. (66) and (81), one can rewrite the definition of the quark condensate (83) in the 
form 

Cquark = 4:N, J ^^cos2vu{q). (84) 

Let us assume that the representation for the wave function Li (81) is still valid for non- 
zero but small quark masses. Then the subtraction of the BS equation (80) from the SD one 
(79) multiplied by the factor 1/-F^ determines the second meson wave function L2 

-Y^I{V) = jr- (85) 

The wave function -^^2 (p) independent of the momentum in this aproximation. Substituting 
the equation L2 = const = 2mu/(M^F^) into the normalization condition (82), and using 
Eqs. (81) and (84), we arrive at the Gell-Mann-Oakes-Renner relation [36] 

M^F^ = 2m,Cq,ark. (86) 

Our solutions including the GMOR relation (86) differ from the accepted ones [30-35], where 
cos2vi^{q) is relaced by the sum of two Goldstone bosons, the pseudoscalar and the scalar 
one [cos2va{q) + (q/o') sin2wa(q')]. This replacement can hardly been justified, because it is 
in contradiction with the Bethe-Salpeter equation (B16) for scalar bound state with nonzero 
mass. 

The coupled equations (71), (74), and (75) contain the Goldstone mode that accompa- 
nies spontaneous breakdown of chiral symmetry. Thus, in the framework of instantaneous 
interaction we prove the Goldstone theorem in the bilocal variant, and the GMOR relation 
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directly results from the existence of the gluon and quark condensates. Strictly speaking, 
the postulate that the finiteness of the gluon and quark condensates are finite implies that 
QCD is the theory without ultraviolet divergence. They can be removed by the Casimir type 
substraction [29] with the finite renormalization [37]. 



6.4. New Hamiltonian interaction inspired by the anomalous triangle diagram with a 

pseudoscalar bound state 

It was shown [22, 23] that the Habbard-Stratanovich linearization of the four fermion 
interaction leads to an effective action for bound states in any gauge theory. We include 
here an effective action describing the direct pion-positronium creation 

where a = 1/137 is the QED coupling constant, and Fp, contained in Eq. (40), plays a 
role of the pion weak coupling parameter = 92 GeV. The first term — ( h ) AiBi 

TT FpJ 

comes from the triangle diagram (i.e., the anomalous term). This term describes the two 7 
decay of pseudoscalar bound states Pbs- 

For each bound state one can obtain the corresponding two-photon anomalous creation 
cross section from Eq. (43). In the case of the process 7 + 7 = Pbs + -Pbs we repeat Eq. (43) 



da a^Et^ / _ / 2m, 



dVL 7rl28F^ V \ 

bs T ^ ' 

where Pbs is the F^^ analogy. In the case of the process 7 + 7 = P^^ + Ppos we obtain 




da a^E^^ 



dn 7r32p2 ^2 



(89) 



The Hamiltonian of this system is the sum of the energy of the free EM fields, the 
pseudoscalar Hamiltonians and their interactions 

W,ff = J dtd\ l^.i, - ?i±^ _ u,^)^ ^ (90) 

This action contains the additional terms in comparison with the standard QED. They leads 
to the additional mass of the pseudoscalar bosons [38] and the anomalous processes of the 
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creation of the bound state pairs in the external magnetic field. The last term of the effective 
action (91) yields cross-sections of creation of both the two positronium atoms and the pion 
and the positronium together. 

The creation of two positronium atoms is times less than the creation of the pion 
and the positronium together. In this case, one can speak about the pion catalysis of the 
positronium creation. 

7. SUMMARY 

In this paper we obtain the bound state functional by Poincare-invariant generalization 
of the FP path integral based on the Markov- Yukawa constraint for description of both the 
spectrum equations and the S-matrix elements. The axiomatic approach to gauge theories 
presented here allows us to construct the bound state functional in both QED and QCD on 
equal footing of the Poincare group representations. 

It is shown that the Poincare S-matrix, as compared with the Lorentz one, contains 

1. Creation of bound states inspired by the anomalous (triangle) diagram within the 
Hamiltonian approach. 

2. This additional anomalous contribution includes the processes like 7 -|- 7 — > Ps + Ps, 
7 -|- 7 — )■ TTo + Ps, 7 -|- 7 — > TTo + ttq (whcrc Ps - a pseudoscalar para-positronium). 

This raises the problem of physical consequences of these additional processes. 

The bound state generating functional (51), where the time-axis is chosen as eigenvalue of 
the total momentum operator of instantaneous bound states (57), has a variety of properties. 
It describes spontaneous breakdown of chiral symmetry, the bilocal variant of the Goldstone 
theorem, and the direct derivation of the GMOR relation directly from the fact of existence 
of the finite gluon and quark condensates introduced by the normal ordering of the QCD 
action. The postulate of the finiteness of the gluon and quark condensates implies that both 
the QED and QCD can be considered on equal footing as the theory without ultraviolet 
divergences. They can be removed by the Casimir type substraction [29] with the finite 
renormalization [37]. 
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APPENDIX A: LADDER APPROXIMATION 



The generating functional (54) can be presented by means of the relativistic generalization 
of the Hubbard- Stratonovich (HS) transformation [22] 



/ + 00 
dyexp[-ixy - y^/i2a)]. 
-oo 



(Al) 



The basic idea of the HS transformation is to reformulate a system of particles interact- 
ing through two-body potentials into a system of independent particles interacting with a 
fluctuating field. It is used to convert a particle theory into its respective field theory by 
linearizing the density operator in the many-body interaction term of the Hamiltonian and 
introducing a scalar auxiliary field [22] 



n dM''\x,y) 

.x,y,a,b 



exp{iWes[M, A*] + i{7]T], Gm)}- 

The effective action in Eq. (A3) can be decomposed in the form 

oo 

W,ff[M, A*] = --N,{M, }C~'M) + iN, -"^^ 

n=l 

Here $ = Ga*-M., etc. mean the following expressions 

$(a;, y) = Ga*M = j d^zGA* {x, z)M{z, y), 
$^ = y d'^xd'^y^{x,y)^{y,x), 
<l>^ = y d'^xd'^yd'^z^{x,y)^{y,z)^{z,x) ,etc 



(A2) 



(A3) 



(A4) 



(A5) 



The first step to the semi-classical quantization of this construction [39] is the determi 
nation of its minimum of the effective action 

_,6W,ff{M) . . i 



n: 



6M 



G^l - M 



0. 



(A6) 
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We denote the corresponding classical solution for the bilocal field by S(x — y). It depends 
only on the difference x — y at A* = because of translation invariance of vacuum solutions. 
The next step is the expansion of the effective action around the point of minimum 

m = j: + m', 

.(2) 



oo oo ^ 

= E ^^"^ = E -iGj:MT, (Ge = {G^l - ^r'), (A7) 



int 

■ ' ■ ' n 

n=3 n=3 



and the representation of the small fluctuations M.' as a sum (5) 

M{x,y) = M{z\X) (A8) 

over the complete set of orthonormalized solutions F, of the classical equation 

^5k£(£±^.r.O (A9) 
with a set of quantum numbers (H) including masses = \/V^ and energies Uh = 



* + M^. The bound state creation and annihilation operators obey the commutation 
relations 



SH'HS'iV - V) . (AlO) 



The corresponding Green function takes the form 

[ [Vo-uh- ie)2ujH [Vo - ujh - t£)2uH J 

To normalize vertex functions, we can use the "free" part of the effective action (A7) for 
the quantum bilocal meson Ai' with the commutation relations (AlO). The substitution of 
the off - shell Vv^ ^ Mh decomposition (5) into the "free" part of effective action deflnes 
the reverse Green function of the bilocal field Q{Vo) 

W^%[M] = 2n6{V, - V'o) J2 J J=«^(^) aJ,{V)g],\Vo) (A12) 

where Qjj^(Vo) is the reverse Green function which can be represented as a difference of two 
terms 

V^\Vo) = I{VV^) - mfiu:)) (A13) 
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where M'^{u) is the eigenvalue of the equation for small fluctuations (A9) and 



X 



tr 



where 



GE,(g-|)ffb( 



1 



^)GE.(g + |)rf,(g^|P) 



is the fermion Green function. The normalization condition is defined by the formula 



(AM) 



2uj 



Finally, we get that solutions of equation (A9) satisfy the normalization condition [40] 



rfM(Po) dI{M) 



\'Po=u} 



(A15) 



d 



d^q 



■tr 



V)G^{q+'^)TH{q^\V) 



2uj 



H- 



(A16) 



dV^j (27r)4 

The achievement of the relativistic covariant constraint-shell quantization of gauge theo- 
ries is the description of both the spectrum of bound states and their S-matrix elements. 

It is convenient to write the relativistic-invariant matrix elements for the action (A7) in 
terms of the field operator 

<D'(a;, y) = j d^xiGj:{x - xi)M'{xi, y) = <^'{z\X) 

Using the decomposition over the bound state quantum numbers [H) 



(A17) 



where 

$H(ab)(g^|P) = Gj,a{q + V/2)THi,i,){q^\V) , (A18) 

we can write the matrix elements for the interaction 14^(") (A7) between the vacuum and the 
n-bound state [20] 



(fflPl...,ff,.P„|!W'"'|0) 



nl/2 



M(")(Pi,...,P„), (A19) 
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{ik} 



where {{ik} denotes permutations over ik). 

Expressions (All), (A17), (A19), and (A20) represent Feynman rules for the construction 
of a quantum field theory with the action (A7) in terms of bilocal fields. 



APPENDIX B: BETHE - SALPETER EQUATION 

Equations for the spectrum of the bound states (A9) can be rewritten in the form of the 
Bethe-Salpeter (BS) one [5] 



r = ilC{x, y) J zidt" Z2Gt.{x - zi)V{zi, Z2)Gt.{z2 - y) . (Bl) 
In the momentum space with 

V{q\V) = j d^xdSe'^^e'^^-y^^Vix, y) 
we obtain the following equation of the vertex function ( F ) 



G^[q+'L)T{q\V)G^{q-'^) 



(B2) 



where V_{k'^) Cleans the Fourier transform of the potential, kj[ = k^ — £^{k ■ i) is the relative 
momentum transversal with respect to i^, and is the total momentum. 
The quantity F depends only on the transversal momentum 

r(fc|P) = r(A:^|p), 

because of the instantaneous form of the potential V_{k-^) in any frame. 

We consider the Bethe - Salpeter equation (Bl) after integration over the longitudinal 
momentum go- The vertex function takes the form 

TAk^lV) = J - 9^)^ab(g^)/, (B3) 
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where the bound state wave function \l/ab is given by 



^ab(g' 



'A 



r,b(g^|P)A(-)b(g^) , A(^)a(g^r,b(g^|P)A(+)b(g 



Et - + it Et + VV^ - ie 

Et = Eg, + E]^ means the sum of one-particle energies of the two particles (a) and 
by (70) and the notation (69) 



A(±)(g^) = 5-i(g^)A(±)(0)5(g^) = A(±)(-g 



(B4) 
defined 

(B5) 



has been introduced. 

Acting with the operators (69) and (B5) on equation (B3) one gets the equations for the 
wave function ip in an arbitrary moving reference frame 



(B6) 



^ti)a(k') I ^Z(^^-g^)v^ab(g-)]Ag),(-^-). 



All these equations (B3) and (B6) have been derived without any assumption about the 
smallness of the relative momentum \k-^\ and for an arbitrary total momentum 



We expand the function \E' on the projection operators 

^ = ^+ + ^_, ^± = Ai^^^A5 
According to Eq. (B4), \E' satisfies the identities 

Ai^^v^Ai^) = aL'^vI/aL') ^ , 



(B7) 



(B8) 



which permit the determination of an unambiguous expansion of in terms of the Lorentz 
structures: 

^a,b± = S;' {75^a,b±(g^) + (7m - A?(0)5b ^ , (B9) 

where L± = Li± L2, N± = Ni± N2. In the rest frame = (1, 0, 0, 0) we get 

= (0, N') ; N\q) = N^{q)eUq) + J:{q)q^ . 

a=l,2 

The wave functions L, N°', S satisfy the following equations. 



1. Pseudoscalar particles. 
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of d?q _ _ _ 

Ml L2 (p) = E Li{v)- j (^^(P ~ q)(S S " ) (^) ! 

(BIO) 

Mi Li (p) = EL2{V)- j (0^^(P - q)(« - ^2 (q) ■ 

Here, in all equations, we use the following definitions 

E(p) = E,(p) + Eb(p) , (Bll) 
c±(p) = cos[vM±Mp)] , (B12) 
s=^(p) = sm[v^{p) ±vi,{p)] , (B13) 
e = ■ , (B14) 

where E^^, E\^ are one-particle energies and v^^, t;b are the Foldy-Wouthuysen angles of particles 
(a,b) given by Eqs.(71) and (72). 

2. Vector particles. 



Mn iV2 " = ^ iVi " - 

/d^n ^ 

^V(p - q){(c;c-r^ + s;s-(r^e - trf)) N, +(r/"s <) Si} ; 

(B15) 

Mn iVi " = ^ iV2 " - 

/^3„ Q /3 Q 

r/" = g,enp) , f = Pi&t{q) , = er(g)ef (p) . 
3. Scalar particles. 




Ms S2 = ^ Si - 



^''^ i^(p - q){(e« + stst) k +{v'c;ct) m 



(2vr) 



(B16) 




Ms El = ^ S2 



"^"^ V{p- q){(ec-c- + S-S-) S2 Hv^ctc-) N2 } . 



(27r) 
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The normalization of these solutions is uniquely determined by equation (A16) 



mJ (2^ 

27V. r d^q 



{LMLl{q) + L,{q)Ll{q)] 



Mn J (27r) 
2Nr f d^q 



3 {iVf (g)iVr (g) + N^{q)Nr{q)] = 1 
' {Si(g)S;(g) + E2(g)Et(g)} = l . 



(B17) 
(B18) 

Ms J (27r)3 ' "-v-^^-iv-^J - • 

If the atom is at rest ( = (Af^, 0, 0, 0) ) equation (B6) coincides with the Salpeter 
equation [4]. If one assumes that the current mass mP is much larger than the relative 
momentum Ig-*"!, then the coupled equations (B3) and (B6) turn into the Schrodinger 
equation. In the rest frame ( Vq = Ma) equation (70) for a large mass {m^/\q-^\ — > oo) 
describes a nonrelativistic particle 

1 k2 



tan2v = ^ ^ 0; ^(k) ^ 1; A(±) ^ 
Then, in equation (B6) only the state with positive energy remains 



^ ~ ^(+) = A(+)75V4/i^/'Sch, A(_)V^A(+) ~ 0, 
where /i = ma ■ rrib/ (rria + mb). And finally the Schrodinger equation results in 



— k" 
2/i 



{ml + ml- Ma) 



fe(k) 



d^q 
(2^ 



Z(k-q)V^Sch(q), 



(B20) 



with the normalization J d^qlipschl"^ /{2 



1. 



For an arbitrary total momentum equation (B20) takes the form 



^iK)~' + iml + m',-Vv^) 



(2vr) 



rH^^-g^)^Sch(g^), (B21) 



and describes a relativistic atom with nonrelativistic relative momentum \k^\ ^ n^ab- 
the framework of such a derivation of the Schrodinger equation it is sufficient to define the 
total coordinate as X = {x + y)/2, independently of the magnitude of the masses of the two 
particles forming an atom. 
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